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CHARACTERISTICS  OF  TURBULENCE  IN  A  BOUNDARY  LAYER 
WITH  ZERO  PRESSURE  GRADIENT 1 

By  I*.  8.  Ki.kbasoff 


SUMMARY 

The  result s  of  an  experimental  investigation  of  a  turbulent 
boundary  layer  with  zero  pressure  gradient  are  presented. 
Measurements  with  the  hot-wire  anemometer  were  made  of 
turbulent  energy  and  turbulent  shear  stress ,  probability  density 
and  flattening  factor  of  u-fluctuation  ( fluctuation  in  x-direction), 
spectra  of  turbulent  energy  and  shear  stress ,  and  turbulent 
dissipation.  The  importance  of  the  region  near  the  wall  and 
the  inadequacy  of  the  concept  of  local  isotropy  are  demonstrated. 
Attention  is  gicen  to  the  energy  balance  and  the  intermittent 
character  of  the  outer  region  of  the  boundary  layer.  Also 
several  interesting  features  of  the  spectral  distribution  of  the 
turbulent  motions  are  discussed. 

INTRODUCTION 

Tin*  statistical  theory  of  turbulence  introduced  by  Taylor 
and  elaborated  upon  by  Kantian,  Howarth,  and  Dryden  has 
played  an  important  role  in  providing  a  sound  basis  for  tlie 
study  of  turbulence.  In  the  main,  this  advance  has  been 
confined  to  the  domain  of  homogeneous  and  isotropic 
turbulence.  At  present  turbulent  shear  flows  present  diffi¬ 
culties  so  formidable  that  statistical  theories  have  made  little 
progress.  The  older  semiempirical  and  phenomenological 
theories  still  constitute  the  most  tangible  theoretical  methods. 

It  is  now  generally  believed  that  experiment  should  be  called 
upon  whenever  possible  to  furnish  information  on  the  actual 
behavior  on  the  grounds  that  such  information  is  needed  to 
acquire  an  insight  into  the  turbulence  processes  and  to  form 
the  basis  for  a  sound  theoretical  approach. 

Experimental  investigations  in  shear  flow  have  made 
encouraging  progress,  largely  because  of  the  increasing  num¬ 
ber  of  statistical  properties  that  may  be  measured  by  im¬ 
proved  hot-wire  techniques.  Measurements  have  been 
made  in  the  jet,  wake,  two-dimensional  channel,  pipe,  and 
boundary  layer  (e.  g.,  refs.  1  to  8).  These  have  shown  that 
the  basic  assumptions  of  the  phenomenological  theories  arc 
inconsistent  with  the  experimental  evidence  as  to  the  nature  i 
of  the  turbulent  motions.  The  work  of  t'orrsin  (ref.  1)  I 
revealed  the  intermittent  character  of  turbulent  How  near  a  ! 
free  boundary.  This  was  later  studied  in  some  detail  by  ! 
Townsend  (ref.  2)  and  is  now  recognized  as  a  phenomenon 
associated  with  a  sharp  but  irregular  division  between  tur¬ 
bulent  attd  nonturbulent  flow.  The  recent  theoretical 
contributions  of  KolmogorofI  (refs.  H  and  10),  Heisenberg 
(ref.  II),  and  others  dealing  with  the  concept  of  local  iso¬ 
tropy  and  the  spectrum  of  turbulent  energy  have  given  added 


impetus  to  experimental  work  in  shear  turbulence  and  have 
encouraged  the  point  of  view  that  some  of  the  properties  of 
isotropic  turbulence  may  be  applicable.  The  application 
of  isotropy  to  shear  flow  has  been  rather  extensively  studied 
in  wakes  (ref.  3)  and  found  to  be  useful,  but  its  usefulness  in 
a  boundary  layer  is  still  questionable,  especially  in  the 
region  near  the  wall.  At  best,  isotropy  belongs  to  the  final 
or  dissipation  stage  of  shear  turbulence  and  here  can  throw 
little  light  on  the  mechanisms  pertaining  to  shear  flow.  Any 
attempt  to  investigate  the  turbulence  mechanism  is  ham¬ 
pered  by  the  lack  of  an  experimental  technique  for  measure¬ 
ment  of  the  pressure  fluctuations.  However,  it  is  possible 
to  obtain  information  pertaining  to  the  energy  balance,  the 
study  of  which  may  be  considered  a  proper  approach  to  the 
problem.  An  important  step  in  this  direction  in  the  investi¬ 
gation  of  the  boundary  layer  was  made  by  Townsend  (ref.  <>). 
With  respect  to  this  point  of  view  the  present  work  attempts 
to  provide  additional  information,  especially  for  the  region 
close  to  the  wall,  and  to  obtain  a  more  direct  measure  of  the 
dissipation  without  relying  completely  on  the  concept  of  local 
isotropy. 

Attention  here  is  given  to  the  boundary  layer  with  zero 
pressure  gradient  under  conditions  as  favorable  ns  practicable 
for  the  use  of  hot-wire  technique.  The  method  of  obtaining 
the  boundary  layer  is  described  in  reference  8.  and  the  present 
work  may  be  regarded  as  a  continuation  of  the  former  work 
in  which  use  is  now  made  of  an  artificially  thickened  layer  as 
a  research  tool.  While  there  are  many  important  features 
of  boundary-layer  turbulence,  (he  present  investigation  is 
concerned  with  three  major  phases:  The  intermit tenev  and 
its  effects  at  the  free  boundary,  energy  balance,  and  the 
spectral  distribution  of  turbulent  energy  and  shear  stress. 
Since  the  boundary  layer  is  a  complex  flow  combining  the 
effects  of  a  fret*  boundary  on  tin*  one  side  and  those  of  a  solid 
wall  on  the  other,  no  one  part  of  the  layer  could  be  do- 
emphasized.  This  makes  it  an  interesting  subject  if  not  a 
simple  one. 

The  present  investigation  was  conducted  at  the  National 
Bureuu  of  Standards  under  the  sponsorship  and  with  the 
financial  assistance  of  tin*  National  Advisory  ('ommillee  for 
Aeronautics.  The  author  wishes  to  express  his  appreciation 
to  l)r,  («.  II.  Sehubauer  for  his  active  support  and  construc¬ 
tive  criticism  and  to  acknowledge  gratefully  the  assistance 
of  Miss  Z.  W.  Diehl  and  Mr.  T.  .1.  Kelly.  He  also  wishes  to 
thank  Dr.  ('.  M.  Tclieti  and  Dr.  .1,  Hauler  for  their  interest 
and  stimulating  discussions. 
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SYMBOLS 

r  45  r"c  ,,  n  ,1  /?+?+p\ 

(’/  coefficient  of  skin  friction,  r„  pi 

DmL\*iy'(u*+lf+*+li) 

Elk)  three-dimensional  spectral  function  associated 

with  k 

e  voltage  fluctuation 

signal  produced  by  H-component 
r.  signal  produced  by  /-component 

F  flattening  factor  _ 

F„(k ,)  percent  of  turbulent  energy  ir  associated  with  k, 

f„lk,)  turbulent  energy  «-  associated  with  k, 

/-(*,)  turbulent  shear  stress  m7  associated  with  /•, 
f,(ki)  turbulent  energy  r  associated  with  kt 

k  three-dimensional  wave  number 

k.  —  ftx3*/***)'1* 

k !  one-dimensional  wave  number 

n  frequency 

l’(u)  probability  density  of  it 

/'( u  it’)  probability  density  of  m/m' 

677?  yr 

’  rt  '  ih, 

/i  mean  static  pressure 

l>„  mean  static  pressure  at  x=  lO.o  feet 

/i,  instantaneous  pressure  fluctuation 

i/„  reference  dynamic  pressure  ahead  of  plate 

Ur  longitudinal  space  correlation  coefficient  of 

it  -fluctuation 

l!„,  lateral  space  correlation  coeflicienls  of  M-fluct na¬ 

tion 

.S'  skewness  factor 

t  time 

I  /-component  of  mean  velocity 


mean  velocity  in  free  stream 

instantaneous  turbulent  velocity  fluctuations 
in  s-, and  redirections,  respectively 
root-mean-squnre  values  of  m.  r,  and  ir 
mean-s(|uare  values  of  m.  r,  and  »• 
turbulent  shearing  stress 
//-component  of  mean  velocity 


...  _  iv  if1  hr  *  r  1  M'*' \ 

2  ) 

(SW+fD'+Q'+O] 

,r  distance  along  surface  from  leading  edge  of  plate 

>/',  intervals  in  /-.  //-,  and  r-direclions 


<listance  along  surface  from  virtual  origin  of 
boundary  layer 

distance  normal  to  surface  measured  from  sur¬ 
face 


direction  perpendicular  to  ./M-plane 
intermit teney  factor 
boundary-layer  t  hickness 
bounds  rv-lavcr  energy  thickness. 


J..  / 


f  dissipation  of  turbulent  energy  in  isotropic 

I  tirbulcncc 
i-=U'2<x  S)- '  |v  5)  rt.TSj 

6  angle  between  velocity  vector  and  normal  to 

wire 

X..  longitudinal  microscale 

v  kinematic  viscosity 

p  density  of  air 

<r  standard  deviation 

t„  shearing  stress  at  wall 

x  absolute  constant 

<j>,  /-component  of  vorticity 

(  )  mean  value 

PROCEDURE  AND  RESl'LTS 
KXl'KMMKNTAI.  ARKANUEMKNT 

The  present  investigation  was  conducted  at  the  National 
Bureau  of  Standards  in  the  4 ‘-.-fool  wind  tunnel  shown  in 
figure  I.  The  turbulence  level  of  the  tunnel  was  (1,02  per¬ 
cent  at  BI I  feet  per  second  and  (1.04  percent  at  Hid  fool  per 
second.  The  low  level  of  turbulence  was  obtained  by  damp¬ 
ing  screens  placed  in  the  settling  chamber.  The  boundary 
layer  was  developed  on  a  smooth,  flat,  aluminum  plate  12 
feet  long.  4'2  feet  wide,  and  h  inch  thick  with  a  symmetrical 
and  pointed  leading  edge.  The  plate  was  mounted  v  ert  ically 
and  centrally  in  the  test  section  of  the  tunnel. 

The  scheme  of  artificially  thickening  of  turbulent  bound¬ 
ary  layers  developed  in  reference  S  was  applied  herein  order 
to  realize  the  advantages  of  a  thick  boundary  layer,  namely, 
the  larger  scales  of  mean  and  fluctuating  How  fields  which 
decrease  errors  due  to  the  finite  size  of  the  hot-wire  probe 
and  the  limited  upper-frequency  response  of  the  equipment 
In  addition,  high  Reynolds  numbers  are  afforded  without 
high  speeds.  The  thickening  was  achieved  by  covering  the 
first  2  feet  of  the  plate  with  sand  roughness  consisting  of 


Sloping  ceiling 
Upper  return  duct 


cone 


Oection  of  flo* 
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No.  16  floor-sanding  paper.  At  the  working  station  IOf* 
feet  from  the  leading  edge  the  boundary  layer  was  :i  inches 
thick.  All  of  the  measurements  were  made  at  this  posi¬ 
tion.  The  free-stream  speed  was  r>0  feet  per  second.  The 
elaborate  tests  described  in  reference  S  showed  that  the 
boundary  layer  was  the  fully  developed,  smooth-wall  type, 
having  an  apparent  development  length  of  14.2  feet  of  smooth 
surface.  The  corresponding  length  Reynolds  number  based 
on  values  of  jr  measured  from  the  virtual  origin  was  4.2X  1 0". 

Ill  order  to  obtain  a  condition  of  zero  pressure  gradient 
along  the  plate,  the  passage  between  the  tunnel  wall  anil  ] 
the  plate  was  made  sufficiently  divergent  to  offset  the  na¬ 
tural  fall  in  pressure  due  to  boundary-layer  growth.  This 
was  accomplished  by  a  flexible  side  wall  which  could  be 
positioned  by  screws  threaded  into  the  tunnel  wall.  The 
final  pressure  distribution  is  shown  in  figure  2.  From  2  feet 
on  downstream  the  pressure  is  seen  to  be  constant  on  the 
average  with  variations  about  %  percent  of  q„  about  the 
mean.  These  variations  were  associated  with  the  inherent 
waviness  of  the  surface  and  could  not  be  removed  by  the 
adjustable  wall. 


The  static-pressure  tube  was  made  according  to  the  con¬ 
ventional  design  for  such  a  tube,  and  the  pressure  distribu¬ 
tion  was  measured  by  traversing  longitudinally  at  u  dis¬ 
tance  *4  inch  from  the  surface.  The  impact  tube  was  flat- 
timed  at  the  end  to  form  a  rectangular  opening  0.014  inch 
wide  and,  together  with  a  static  tube  similar  to  that  used 
for  the  pressure  distribution,  was  used  to  measure  the  mean- 
velocity  distribution.  Velocity  distributions  were  also 
measured  at  10  inches  above  and  below  the  center  line. 
These  agreed  well  with  the  distribution  at  the  center  line 
and  thus  confirmed  the  two-dimensional  nature  of  the  flow. 
No  correction  was  made  for  the  effect  of  turbulence  on  the 
measured  values  nor  for  the  effective  geometric  center  of  the 
impact  tube. 

(’lose  to  the  wall.  0<  i/^  II. On  inch,  the  mean  velocity  was 
measured  with  a  hot-wire.  Platinum  wire  O.OOO:!  inch  in 
diameter  and  approximately  L  inch  long  was  used  and 
operated  at  low-temperature  loadings  in  order  to  minimize 
tin1  influence  of  the  wall  on  the  heat-loss  characteristics  of 
the  wire.  A  correction  was  made  for  the  effect  id'  the  turbu¬ 
lence  level  on  the  measured  hot-wire  values.  This  was  done 


INSTRUMENT  MOUNTING  AND  TRAVERSING 

The  positioning  and  moving  of  an  instrument  from  point 
to  point  was  always  done  in  the  manner  that  best  suited  the 
instrument  and  the  purpose.  For  example,  the  pressure 
distribution  along  the  plate  was  determined  by  means  of 
a  static  tube  mounted  from  a  carriage  that  could  be  moved 
and  positioned  by  remote  control.  Since  all  other  types  of 
measurements  were  made  at  the  lO.n-foot  station,  the  various 
measuring  probes  were  supported  on  rods  extending  through 
the  plate  to  a  micrometer-screw  traversing  device  on  the 
opposite  side.  This  provided  adequate  rigidity  and  negli¬ 
gible  interference  and  permitted  movement  by  known 
amounts  to  and  from  the  surface.  Initial  distance  from  the 
surface  was  obtained  by  using  a  prism  to  reflect  the  images 
of  the  surface  and  the  probe  on  the  calibrated  scale  of  a 
microscope. 

In  special  eases  such  as  the  measurements  of  the  //,  and 
ll:  correlations,  small  manually  operated  traversing  units 
that  could  be  mounted  to  the  rods  were  employed.  In  all 
eases  adequate  rigidity  and  freedom  from  interference  rather 
than  convenience  of  operation  dictated  the  arrangement. 

MEASUREMENT  OF  PRESSURE  AND  MEAN  VELOCITY 

Nickel  t tilling  0.04  inch  in  diameter  and  O.OOtt  inch  in  wall 
thickness  was  used  for  the  impact  and  static-pressure  tubes. 
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by  a  graphical  method  using  the  mean-veloeit  \  -volt age  cali¬ 
bration  curve  and  the  measured  root  mean  square  of  the 
voltage  fluctuations.  The  corrected  mean  velocity  was 
higher  than  the  observed,  with  a  maximum  correction  of 
about  lit  percent.  The  hot-wire  values  were  in  good  agree¬ 
ment  with  the  pilot  -st  a  t  ie-t  it  he  values. 

'Pile  velocity  distribution  is  shown  in  figure  The  dashed 
line  denoting  the  velocity  gradient  at  I  lie  wall  was  computed 
from  the  shear  at  the  wall  and  is  in  satisfactory  agreement 
with  t he  measured  values. 

TURBULENCE  INTENSITIES  AND  SHEARING  STRESS 

The  hot-wire  equipment  used  is  described  in  detail  in 
reference  12.  The  frequency  response  of  the  uncompensated 
amplifier  was  flat  from  2  to  of), 01  Ml  cycles  and  was  down  17 
percent  at  70.000  cycles.  number  of  cutoff  filters  were 
provided  to  cut  out  frequencies  above  the  range  needed  in  a 
particular  measurement  so  that  unnecessary  noise  was 
eliminated.  By  t his  expedient  the  input  noise  level  was  held 
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down  to  the  order  of  2  to  4  microvolts.  Compensation  for 
the  time  lag  of  the  wire  was  determined  by  the  square-wave 
method  and  was  accomplished  by  a  resistance-capacitance 
network  in  the  amplifier.  Platinum  and  platinum-rhodium 
(90  percent  platinum  and  10  percent  rhodium)  wires  0.0001 
inch  in  diameter  and  %  millimeter  long  having  a  time  constant 
of  approximately  0.25  millisecond  were  generally  used.  In 
cases  where  the  signal-to-noise  level  presented  some  difficulty, 
0.00005-inch-diameter  wire  was  used.  This  was  wire  drawn 
by  the  Wollaston  process  and  the  method  of  attaching  the 
wire  was  to  etch  away  the  silver  and  to  soft-solder  the  wire 
to  the  supporting  prongs.  Fine  sewing  needles  were  used  as 
the  prongs  for  the  u-holder  anti  fine  jeweler’s  broaches  were 
used  as  prongs  for  the  X-wires.  Xo  wire-length  corrections 
were  applied  to  any  of  the  data.  The  method  of  measuring 
the  turbulence  intensities  r\  and  w'  is  described  in  refer¬ 
ence  Id.  The  i/’.  r',  and  i r'  distributions  are  shown  in  figure 
4.  The  value  of  u'  was  obtained  as  close  to  the  wall  as  0.004 
inch,  but  because  of  the  comparatively  larger  size  of  the 
probes  necessary  for  the  measurement  of  r'  and  «*'  it  was 
not  possible  to  measure  these  closer  than  about  0.045  inch. 
The  extrapolated  values  of  r'  and  w'  were  obtained  by  com¬ 
paring  the  boundary-layer  values  with  those  obtained  in  a 
pipe  (ref.  5)  on  a  basis  of  l  \  against  y*.  Good  agreement 
existed  within  the  measured  range  and  consequently  the  pipe 
data  which  contained  values  corresponding  to  distances  closer 
to  the  wall  (because  of  the  lower  Reynolds  number)  were  used 
as  a  basis  for  extrapolation. 

The  measurement  of  577  was  made  by  the  conventional 
X-wire  method  described  in  reference  III.  Cse  was  made  of 
the  experimental  data  for  the  angle  response  found  by 
Newman  and  Leary  (ref.  14).  namely,  (cos  0)”  ,5;  rather  than 
the  customary  (cos  6)],i  where  8  is  tin*  angle  between  the 
velocity  vector  and  the  normal  to  the  wire.  The  signals  to 
be  dealt  with  are  fx  and  e2  produced  by  the  u-  and  /--com¬ 
ponents,  respectively.  In  principle,  one  wire  of  an  X-wire 
probe  contributes  the  output  (c,i-r2)‘  ami  the  other  con¬ 
tributes  (e |  —  e2)‘.  The  difference  between  these  two  outputs 


gives  7j72  from  which  377  may  be  calculated  by  employing 
the  experimentally  determined  factors  of  proportionality. 

The  results  are  shown  in  figure  5.  Since  the  viscous  shear¬ 
ing  stress  reached  only  2  percent  of  the  total  at  the  point 
nearest  the  wall,  the  turbulent  shearing  stress  377  shows  the 
characteristics  of  the  total.  It  approaches  the  wall  with 
zero  slope,  us  it  must  when  the  pressure  gradient  is  zero  and 
is  in  excellent  agreement  with  the  value  at  the  wall  calculated 
by  the  method  of  Squire  and  Young  (ref.  15). 

SPECTRA  OK  E.  AND  u* 

The  mean-square  values  of  the  u-fluctuation  between  fre¬ 
quencies  n  and  n+ihi  were  obtained  by  feeding  the  signal 
from  the  compensated  amplifier  into  a  General  Radio  wave 
analyzer  having  a  frequency  range  from  10  to  10.000  eps  and 
a  fixed  band  pass.  The  signal  from  the  analyzer  was  then 
fed  into  a  thermocouple  circuit  for  measuring  the  mean- 
square  output.  In  reading  the  output,  averages  were  taken 
j  for  a  period  of  1  minute.  The  menu-square  voltage  asso¬ 
ciated  with  the  hot-wire  signal  passed  by  the  hand  wits  ob- 
i  tabled  by  feeding  in  a  known  sine-wave  input  at  the  proper 
frequency  to  give  the  same  output  reading.  The  mean- 
square  voltage  per  unit  frequency  was  then  obtained  b\ 
dividing  by  the  effective  band  width.  The  effective  band 
width  had  a  value  of  5. .‘10  eps  and  was  defined  as  the  rec¬ 
tangular  band  width  having  the  same  area  as  the  experi¬ 
mentally  determined  band  shape.  The  results  presented  in 
figures  0  and  7  are  normalized  and 


and  the  wave  number  is  given  by 


Over  most  of  the  range  the  accuracy  was  on  the  order  of 
±10  percent.  The  acctiracx  is  somewhat  less  at  the  two 
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riM'KK  H. —  >S|i«»rtra  of  t iirt>ult*iit  shear  stress. 


The  spectra  of  /-’shown  in  figures  1)  anil  1(1  were  calculated 
from  the  data  taken  during  the  measurements  of  the  shear 
spectra.  This  method,  in  principle,  is  perhaps  not  so  accu¬ 
rate  as  a  direct  measurement  of  the  /'-spectra,  hut  a  cheek  on 
the  reliability  of  the  measurement  was  afforded  by  comparing 
the  //-spectra  calculated  from  the  same  data  with  that  di¬ 
rectly  measured,  and  in  general  the  agreement  was  good. 

The  data  for  the  shear  and  r‘  spectra  are  not  normalized 
so  that 


/'X 


f  u>  f^*l )  '/I'l 
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PROBABILITY  OISTRIBI'TION.  SKKWNKSS.  ANI>  FI.ATTKNINU  FACTOR  or 
u-FI.IT  TI  ATIONS 

The  distribution  of  amplitude  of  //-flucl nations  was  de¬ 
termined  by  passing  the  voltage  signal  through  a  eon veu 
lional  gate  circuit  which  conducted  when  the  signal  was 
above  the  level  of  the  gate.  The  signal  was  preamplified  so 
that  the  gate  width  could  be  considered  neglible  with  respect 
to  the  root-mean-scpiare  value  of  the  signal  and  tin-  pulses 
produced  by  the  gate  were  a  series  of  square  waves  of  con¬ 
stant  amplitude  of  varying  duration  dependent  on  the  time 
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Kn.riu;  10  ('omparisun  <if  measurctl  r2  s|H*rtriiin  with  that  (-nictitated  from  «2  spectrum  using  isotropic  relation;  y  &  0  o Y 


varied  with  the  level  of  the  pile  and  ranged  from  I  minute 
at  the  oil-percent  position  to  it  minutes  at  the  edge  where  the 
pulsi>s  were  short  and  infretpient. 

The  midpoint  was  defined  as  the  position  about  which  the 
average  signal  was  zero  and  is  obtained  from  the  lirsl  moment 
of  the  |>robabi)ity  distribution,  which  by  definition  is 

a /’( it )  i/ii  (I 

The  distributions  of  the  probability  density  at  various 
positions  across  the  boundary  Inver  are  shown  in  figure  1 1 
and  are  expressed  in  values  of  u  relative  to  so  that 


J7CMJ) 


At  most  of  (he  positions  the  distributions  closed  in  rather 
well  anti  permitted  the  calculation  with  reasonable  accuracy 
of  the  skewness  and  flattening  factors,  when* 
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At  ;/  4  O.M  ami  y:i-- 1 .0  the  limits  of  the  curves  were  uncer- 
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tain,  and  no  attempt  was  made  to  calculate  the  higher 
moments  at  these  positions.  The  flattening  factor  was  also 
measured  directly.  Since  this  involves  the  nondimcnsional 
ratio  of  the  same  voltage  signal,  that  is, 

F=  =  M7 
(T5)2  (a5)2 

it  was  not  necessary  to  calibrate  the  hot-wire,  but  the  cor¬ 
rection  for  the  difference  in  gain  involved  in  the  processing 
of  the  divided  signals  had  to  be  determined  by  calibrating 
with  a  known  sine-wave  input.  The  fourth  power  was 
measured  by  feeding  the  signal  from  the  compensated  ampli¬ 
fier  into  an  instantaneous  squaring  circuit,  which  consisted 
of  a  series  of  diodes  properly  biased  to  give  the  square  of  the 
input  (ref.  12),  and  then  squaring  again  with  a  thermocouple. 
The  distribution  of  flattening  factors  across  the  boundary 
layer  measured  in  this  manner  is  given  in  figure  12. 


JL 

u' 


Kiiiurk  II.  -Distributions  of  probability  density. 
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PiriVRK  12.  -Distribution  of  flattetiinK  factor. 


It  is  seen  from  figure  11  that  from  »//5=0.4  to  the  wall  tin- 
usual  safety  factor  of  three  times  the  root-mcan-squarc  value 
for  determining  the  operating  point  of  the  equipment  is  quite 
satisfactory.  In  the  outer  region  of  the  layer,  at  ;/  a  (l.s. 
this  factor  is  exceeded  in  the  direction  of  negative  values  of  u . 
Here  measurements  were  taken  at  varying  gain-coutml 
settings  to  make  sure  that  the  equipment  was  not 
overloading. 

MEASUREMENT  OF  DERIVATIVES 

The  usefulness  of  the  concept  of  local  isotropy  lies  in  the 
hope  that  the  rate  of  turbulent-energy  dissipation  in  shear 
flow  can  be  given  in  terms  of  the  dissipation  expression  for 
isotropic  turbulence.  This  would  simplify  the  experimental 
procedure  by  requiring  only  the  measurement  of  the  mean- 
square  derivative  of  the  u-fluctuation  with  respect  to 
A  fuller  discussion  as  to  the  adaptability  of  the  concept  of 
local  isotropy  is  given  in  the  section  entitled  “Knergy 
Balance,”  but.  it  may  be  stated  that  the  use  of  the  isotropic 
relation  in  determining  the  dissipation  is  inadequate, 
especially  in  the  region  near  the  wall.  It  thus  becomes 
necessary  to  measure  all  of  the  following  nine  mean-square 
derivativesappearing  in  the  dissipation  term  of  the  turbulence 
energy  equation 

o-o-o-o-tri- 

(IN!^<7<7 

In  the  present  investigation,  because  of  the  practical 
limitations  of  hot-wire  techniques  it  was  possible  to  measure 
only  tin-  first  five  of  the  above  terms.  The  first  three  were 
evaluated  by  taking  the  time  derivative  of  tin-  signal  b\ 
means  of  a  resistance-capacitance  network  and  obtaining  the 
mean-square  value  with  a  thermocouple  and  meter.  The 
practical  use  of  a  differentiating  circuit  involves  a  compromise 
between  attenuation  and  the  extent  of  linear  response  with 
frequency.  In  the  present  measurements  the  differentiation 
was  linear  to  10,000  eps  and  was  down  12  percent  from 
linearity  at  10,000  eps.  This  was  considered  adequate  for 
the  frequencies  under  investigation  since  no  significant  change 
in  the  measured  values  was  obtained  by  varying  the  cut-off 
filters  in  the  amplifier.  The  time  derivative  was  converted 
to  the  space  derivative  by  assuming  the  accuracy  of  the 
space-time  t ransforma t ion 


The  remaining  two  terms  were  obtained  by  measuring  the 
correlation  coefficients  ]{„  and  It,  for  n  small  distance  of  >/ 
and  ;■  and  using  the  relations  given  by  Taylor  (ref.  17) 


n  M(j/)«(2/+.V')  .  _  1 

u'(y)u'(y  +  y')  ->m2 


and 


J/’ — *0 
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„  _  u(z)u(z+z')  _  1  /  duV  (  , 

*  tt'(2jM'(2+2')_  ajpVdr^''' 


The  correlation  coefficients  were  measured  by  the  method 
described  in  reference  18  for  varying  distances  of  y'  and 
ranging  from  0.004  to  0.04  inch. 

The  distributions  of  the  various  derivatives  across  the 
boundary  layer  are  given  in  figures  13  and  14. 

An  attempt  to  assess  the  validity  of  the  space-time  trans¬ 
formation  was  made  at  y/i=0.05.  In  figure  15,  the  longi¬ 
tudinal  correlation  coefficient  Hx  obtained  by  measuring  the 
correlation  between  values  of  u  at  the  same  instant  for 
varying  distances  in  the  /-direction  is  compared  with  that 
calculated  from  the  spectrum  using  the  Fourier  transform 


"■  J> 


(t |)  COS  ilt-fj-')  (/<•, 


The  latter  gives  the  spatial  correlation  from  a  time  correla¬ 
tion  using  the  space-time  relat ion.  The  microscale  X,  calcu¬ 
lated  from  the  spectrum  is  defined  by 


It  is  seen  that  for  values  of /'  less  than  0.2  inch  the  two  are 
in  good  agreement.  They  yield  the  same  microscale,  and  the 
space-time  transformation  is  apparently  valid  for  the  small- 
scale  motions  responsible  for  the  dissipation.  For  large 
values  of  /'  the  two  begin  to  diverge  and  the  space-time 
relation  becomes  progressively  worse  for  the  larger  scale 
motions.  At  r'~0. 4  inch  the  major  contribution  to  the 
correlation  comes  from  those  wave  numbers  below  I  H  per 
centimeter.  It  seems  therefore  that  the  adequacy  of  the 
space-time  transformation  depends  on  wave  number  and 
gives  rise  to  the  interesting  speculation  that  the  large-scale 
motions  have  their  own  characteristic  velocities  different 
from  the  mean  speed.  Still  closer  to  the  wall,  the  space-time 
assumption  does  not  become  a  serious  obstacle  to  the  calcula- 


°  Measured 

a  Calculated  from  spectrum 


Fiourk  15.  Comparison  of  loiiKitudinal  correlation  coefficient  Ht 
obtained  directly  with  that  eaten  fatted  from  spectrum  of  y/4=*O.05. 
Curve  for  X,  also  calculated  from  spectrum. 
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lion  of  llie  dissipation  since  the  magnitude  is  given  mainly  by 
derivatives  in  the  y-  and  ’-directions. 

Although  no  wire-length  corrections  were  applied  to  any  of 
the  data,  it  should  be  mentioned  that  an  estimate  of  the 

bu 
by 

given  in  reference  19.  It  was  approximately  10  percent  and 
consequently  wire-length  corrections  may  be  considered 
negligible  across  most  of  the  boundary  layer. 

DISCUSSION 

INTERMITTENCY  • 

In  shear  Hows  that  have  a  free  boundary,  it  has  been  re¬ 
peatedly  observed  that  as  the  free  stream  is  approached  the 
turbulence  becomes  intermittent,  that  is,  that  for  only  a 
fraction  Y  of  the  time  is  the  flow  turbulent.  This  on-and-off 
character  of  the  turbulence  has  been  definitely  established  as 
being  a  manifestation  of  the  irregular  outline  of  the  boundary 
layer  as  it  moves  downstream.  The  intermittency  is  easily 
observed  by  oscilloscope  records  of  the  M-fluetuution  in  the 
outer  region  of  the  boundary  layer,  and  the  records  can  be 
used  to  give  a  quantitative  estimate  of  the  factor  y  and  to 
discern  some  qualitative  aspects  of  the  flow.  Representative 
sections  of  oscilloscope  records  taken  at  various  positions 
across  the  boundary  layer  are  given  in  figure  Hi.  It  is  seen 
that  in  the  outer  region  of  tin1  layer  //  5>0.4  there  are  inter¬ 
vals  of  time  when  the  flow  is  not  turbulent  and  that  this  time 
increases  with  increasing  distance  from  the  wall.  Thus,  the 
outer  regime  is  divided  into  a  turbulent  part  and  a  relatively 
nontiirbulcnt  free-stream  part,  and  the  hot-wire  at  a  given 
position  responds  to  alternate  turbulent  and  nonturbulcnt 
Hows  as  the  pattern  is  swept  downstream.  A  vortieitv  meter 
i ref.  20)  which  responds  to  the  vortieitv  that  exists  only  in 
the  turbulent  region  and  suppresses  the  low-frequency 
fluctuation  was  used  to  obtain  a  sharper  division  between  the 
turbulent  nnd  nonturbulcnt  regions.  Typical  samples  of  the 
vortieitv  trace  uz  at  different  positions  across  the  boundary 
layer  arc  shown  in  figure  17. 

*  MriiSUirmi-nts  by  Corrsin  ami  Kisth*r  <rrf.  21*  show  that  the*  intormlttunt  region  in 
a  boundary  layer  o\rr  a  rough  wall  has  rhu  sum**  general  behavior. 


KuiCRK  It).  Turbulence  velocity  u  in  boundary  layer.  f'h  50  feet 
per  second;  timing  dot*,  60  per  second. 


correction  was  made  for 


at  y! 6=0.005  by  the  method 


The  first  measurements  of  the  intermittency  factor  y  were 
those  obtained  by  Townsend  (ref.  if),  where  y  is  taken  as  the 
ratio  of  the  flattening  factor  (of  u  or  bulbt)  in  the  wholly 
turbulent  region  to  that  in  the  intermittent  region.  From 
figure  12,  it  is  seen  that  the  flattening  factors  in  the  region 
near  the  wall  y/5<0.4  are  quite  constant  with  a  value  cor¬ 
responding  closely  to  the  Gaussian  value  of  :{.().  By  con¬ 
sidering  the  intermittency  as  an  on-off  process  the  value  of 
y  is  given  by 


The  intermit teneies  calculated  in  this  manner  together  with 
those  calculated  from  the  oscilloscope  records  of  the  //-fluc¬ 
tuation  and  the  vortieitv  are  given  in  figure  lk.  Values  of 
y  obtained  from  the  flattening  factor  are  consistent  with 
those  obtained  from  the  film  only  up  to  y  5-0.9.  As  seen 
in  figure  12,  tbo  flattening  factors  reach  a  maximum  and  then 
begin  to  decrease.  This  is  not  too  surprising  because  for 
low  values  of  7  the  turbulence  of  the  free  stream  would  be 
expected  to  make  itself  felt.  The  result  is  probably  a  weight¬ 
ing  of  the  probability  density  for  the  turbulence  of  the  free 
stream  with  that  within  the  boundary,  which  depends  on 
the  relative  turbulence  levels  and  the  degree  of  intermit tenev. 
The  curve  in  figure  18.  which  closely  represents  the  variation 
of  7  with  yIS,  is  a  Gaussian  integral  curve  given  by 

7=4  (1  ~<'rf  f) 

where 

[(?)-» »] 

and 


The  standard  deviation  a  is 


<r  =  0.145 

Such  a  distribution  indicates  that  the  instantaneous  posi¬ 
tion  of  the  edge  of  the  boundary  layer  lias  a  random  character 


Figure  17. — Turbulence  velocity  u,  in  Iwtuidary  layer.  (  *,  50  feet 
per  second;  liming  dots,  (10  per  second. 
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with  u  mean  position  at  y/S—0.78.  'I'he  edge  rarely  extends 
outside  the  region  y/S= 0.4  to  1.2.  While  the  position  of 
the  edge  fluctuates  over  a  large  fraction  of  6,  the  mean  veloc¬ 
ity  is  near  that  of  the  free  stream,  being  down  by  at  most 
al>out  15  percent. 

Several  interesting  features  can  he  gleaned  from  careful 
study  of  the  oscilloscope  records.  Although  it  may  he  diffi¬ 
cult  to  see  from  the  short  sections  given  in  figure  16,  it  was 
noticed  that  in  the  strongly  intermittent  region,  for  example 
at  y/6= 0.8,  the  trace  had  somewhat  of  a  square-wave  appear¬ 
ance.  Increasing  velocities  are  in  the  direction  of  the  timing 
signal.  The  nonturbulent  regions  seem  to  he  at  a  constant 
level  corresponding  to  that  of  the  free  stream,  while  the 
turbulent  regions  are  seen  to  he  centered  about  some  lower 
level.  The  difference  between  the  velocity  of  the  outside 
potential  flow  and  that  existing  in  the  turbulent  region 
seemed  to  depend  on  how  far  past  the  measuring  position 
the  instantaneous  edge  of  the  layer  extended  at  the  particular 
instant.  At  y/ 5=1.0  and  1.2  there  is  very  little  evidence  of 
the  shift  because  the  edge  does  not  extend  to  any  great  dis¬ 
tance  beyond  these  positions.  There  is  then  a  step  as  well 
as  the  on-off  process  which  was  assumed  in  the  calculation 
of  y  from  the  flattening  factor.  The  agreement  with  the 
values  obtained  from  the  film  is  probably  due  to  the  lack  of 
sensitivity  of  the  flattening  factor  to  this  shift.  However, 
the  skewness  factor  can  be  expected  to  be  extremely  sensitive. 

A  sketch  of  the  boundary  layer  is  given  in  figure  19.  The 
boundary  between  the  turbulent  flow  and  the  free  stream  is 
quite  sharp,  and  the  properties  of  the  shear  layer  an*  com¬ 
paratively  distinct  from  those  in  tin*  free  stream.  The 
boundary  layer  travels  downstream  with  an  outline  con¬ 
stantly  changing  in  an  irregular  manner,  and  the  intermit- 
teney  is  characterized  by  a  large-scale  diffusion  process, 
carrying  with  it  small-scale  turbulent  motions.  From  the 
film  taken  at  y  16=0.8  it  was  also  noticed  that  the  average 
frequency  of  the  occurrence  of  periods  of  no  turbulence' 
seemed  to  be  approximately  100  per  second,  and.  since  tin1 


Fkiuhk  18, — Distribution  of  interinittency  factor. 


13 

pattern  is  moving  with  nearly  the  free-stream  speed  of  50 
feet  per  second,  a  rough  estimate  of  the  average  wave  length 
for  the  irregular  outline  of  the  boundary  layer  would  appear 
to  be  approximately  26  or  1 4<r. 

'I'he  effect  of  the  intermit tency  on  the  probability  density 
is  clearly  seen  from  figure  11.  In  the  noniutermit tent 
region  the  distributions  are  very  nearly  Gaussian  and  the 
values  of  skewness  and  flattening  factor  calculated  from  the 
measured  probabilities  are  given  in  the  following  table: 


o.noi  2  .ra  «.<«> 

.«:»  2.7  f>  —ok 

2  2.U2  — .  OK 

4  j  3.  1H  — .  2»i 

Calculations  of  the  higher  moments  of  the  probability 
densities  from  the  distributions  tend  to  In'  inaccurate  because 
of  their  emphasis  on  the  higher  values  of  it.  However,  the 
;  values  of  flattening  factor  agree  fairly  well  with  the  directly 
measured  values.  Values  of  skewness  should  be  taken  with 
some  reservation  because  they  are  extremely  sensitive  to  the 
accuracy  of  the  midpoint.  Also  in  the  region  near  the  wall 
where  turbulent  fluctuations  are  large  the  nonlinear  response 
of  the  hot-wire  tends  to  skew  the  signal.  At  y  6  O.X  the 
probability  density  is  very  strongly  negatively  skewed 
because  of  the  lower  velocities  within  the  turbulent  regions. 
At  the  same  time  the  large  percent  of  time  that  the  nontnr- 
bnlent  regions  exist  causes  the  maximum  to  be  displaced  to 
the  positive  side  of  the  midpoint.  The  distribution  at  y  6  -  1.(1 
hns  a  flattening  in  qualitative  agreement  with  the  trend  of  the 
flattening  factor  to  reach  a  maximum  and  then  to  decrease. 
The  distributions  of  skewness  show  the  same  general  trend. 
The  maximum  of  the  probability  distribution  for  y  6  1.(1  is 

displaced  slightly  to  the  left  of  »/=-().  and  it  is  uncertain 
whether  this  is  experimental  or  a  consequence  of  the  weight  ing 
effect  previously  mentioned  in  connection  with  the  flattening 
factor. 

By  assuming  that  the  free-stream  regions  contribute  little 
to  the  measured  mcnn-turbuIcncc  quantities,  an  allowance 
may  be  made  for  the  effect  of  interinittency  by  dividing  by 
the  factor  7.  The  distributions  of  turbulent  energy  and 
shear  stress  divided  by  7  are  given  in  figure  20.  The  dis¬ 
tribution  of  turbulent  energy  within  the  bounded  flow  is 
,  strikingly  similar  to  that  for  channel  and  pipe  flow  (refs.  4 
and  5)  as  is  the  distribution  of  turbulent  shear  stress  which 
is  approximately  linear  from  y  6  —  0.1  to  y  d  1.0. 

Interinittency  is,  of  course,  absent  from  fully  developed 


boundary  between  loye; 


Fk.vhk  Ml.  Sketch  of  I  oundarv  layer. 
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turbulent  flow  in  pipes  aiul  channels  because  there  is  no  free 
stream.  It  is  well-known  that  the  mean  velocity  for  pipe 
ami  channel  does  not  deviate  as  much  from  the  logarithmic 
distribution  as  does  that  of  the  boundary  layer,  and  a  major 
factor  in  this  discrepancy  is  the  intermit tencv.  The  in¬ 
fluence  of  the  intermittency  on  the  mean  velocity  distribu¬ 
tion  is  difficult  to  ascertain,  being  in  the  nature  of  a  com¬ 
plex  time-averaging  problem.  Although  something  may  be 
done  in  a  semiempirieal  fashion  by  introducing  7  into  the 
mixing-length  theories,  it  contributes  little  to  the  bsi«c 
understanding  of  the  problem.  In  figure  21  the  measured 
mean  velocity  is  plotted  in  the  form  suggested  by  the  veloc¬ 
ity-defect  law,  and  it  is  seen  that  the  characteristic  loga¬ 
rithmic  law  exists  only  in  a  limited  range  from  about  y/&— 
0.01  to  i//4  =  0.2.  The  direct  viscous  dissipation  (fig.  22)  is 
negligible  at  »//5=0.0l.  This  corresponds  to  a  value  of  y*= 
27.0,  which  is  in  good  agreement  with  the  usually  observed 
value  of  :10  for  pipe  and  channel.  An  interesting  point  to 
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Ficvre  20.  —Turbulent  shear  stress  and  turbulent  energy  divided  by  -y. 


note  is  that  in  the  range  of  j//6=0.1  lo  y;6= O.N  the  ratio  of 
shear  stress  to  turbulent  energy  is  approximately  constant, 
as  is  also  evidenced  by  the  constancy  of  the  shear  correla¬ 
tion  coefficient  in  figure  5.  This  would  tend  to  indicate 
some  degree  of  “Karman  similarity”  existing  in  the  range 
where  the  distribution  deviates  from  the  logarithmic  law. 
In  addition,  K&rmlin  similarity  requires  that  the  large-scale 
motions  responsible  for  the  shear  stress  be  free  from  the  ef¬ 
fects  of  viscosity,  which,  as  will  be  seen  in  the  section  on 
spectra,  is  much  more  justified  in  the  region  where  the 
logarithmic  law  is  not  obeyed  than  in  the  region  win-re  it  is 
obeyed.  Although  any  correction  for  the  effect  of  intermit¬ 
tency  is  in  the  direction  of  minimizing  this  deviation,  it  is 
pointed  out  that  the  mean-velocity  distributions  for  pipe 
and  channel  also  deviate  in  this  region.  It  is  difficult  to 
attach  any  particular  significance  to  the  logarithmic  law 
which  exists  for  so  limited  a  range  of  the  boundary  layer, 
except  that  it  is  a  region  where  the  direct  influence  of  the 
wall  may  slill  be  present.  It  is  not  apparent  that  a  degree 
of  validity  can  be  assigned  lo  it  on  the  basis  of  the  various 
forms  of  the  mixing-length  theory.  With  respect  to  the 
mean  velocity,  the  boundary  layer  can  be  divided  into  three 
!  regions,  namely,  a  viscous  region  extending  to  a  value  of 
y*=:W,  an  intermediate  region  where  the  influence  of  the 
i  wall  still  exists,  ami  an  outer  region  characterized  by 
intermit  tencv. 

I  * 


Flul'RK  21. — Comparison  of  measured  mean  velocity  distribution  with 
velocity-defeet  law. 


1  Fiuckk  22.  Comparison  of  direct  viscous  dissipation  with  production 
of  turbulent  energy  near  wall. 


CHARACTERISTICS  OF  TURBULENCE  IN  A  BOUNDARY  LAYER  WITH  ZERO  PRESSURE  GRADIENT 


lf> 


ENERGY  BALANCE 

After  the  usual  boundary-layer  approximations  are  made, 
the  equation  expressing  the  energy  balance  for  the  turbulence 
at  a  given  cross  section  in  a  two-dimensional  boundury  layer 
is  given  by 

~  %  +2  \  §y  Cw)+l  r  ~  (P+7+^F)  + 

*  V’ ~  (uiJrei+ur)  —  v(uV!u+i'Y2rJr  wV!w)=0 

The  respective  terms  from  left  to  right  have  the  following 
physical  interpretation : 

(1)  Production  of  turbulent  energy  from  the  mean  motion 

(2)  Turbulent  energy  diffusion 

(.'{)  Pressure  diffusion 

(4)  Convection  of  turbulent  energy  by  the  s-componcnt  of 
the  mean  motion 

(5)  Convection  of  turbulent  energy  by  the  (/-component  of 
the  mean  motion 

Hi)  Dissipation  of  turbulent  energy 

As  discussed  by  Corrsin  in  reference  22  the  viscous  term 
in  the  energy-balance  equation  is  not  in  a  strict  sense  the 
true  dissipation  term.  However,  since  for  most  of  the 
boundary  layer  except  for  the  region  close  to  the  wall,  as  well 
as  in  the  integral  across  the  boundary  layer,  it  reflects  the 
true  dissipation,  it  is  herein  referred  to  as  the  dissipation 
term. 

The  convection  of  turbulent  energy  is  negligible  near  the 
wall  and  is  significant  only  in  the  outer  region  of  the  layer 
where  similarity  on  tin1  scale  of  6  has  been  experimentally 
observed.  Consequently,  by  using  the  continuity  relation 

or ,  dl'  n 

dr  d(/ 


aocj-,i45 

the  convection  terms  can  be  transformed  into  the  single  term 
which  may  be  written  in  the  nondimensionnl  form 

_J4  &  fw‘r. 1  <1  fi?+?+P\ 

U  f  ,r„  J„  rrt  ,1{w6)}<k„!6)\  2  ) 

where  i„  is  the  distance  from  the  virtual  origin  of  the  bound¬ 
ary  layer.  The  dissipation  term  is  not  in  a  form  that  lends 
itself  to  measurement  by  hot-wire  techniques.  However,  it 
can  be  rewritten  in  the  following  more  suitable  form: 

All  of  the  pertinent  terms  in  the  energy  balance  can  be 
calculated  from  the  measurements  previously  discussed  ex¬ 


cept  t lit*  diffusion  terms  which  arc  treated  as  one  and  obtained 
by  balancing  the  equation.  The  distributions  of  the  various 
energy  terms  are  given  in  figures  22,  24,  and  25  in  uondimen- 
sional  form  using  5  and  The  turbulent  shear  stress  in 
the  region  ()<j//6< 0.05  was  obtained  by  subtracting  the 
viscous  shear  stress  (calculated  from  the  mean-velocity 
gradient)  from  the  constant  value  of  total  shear  stress  as 
given  by  the  shear  at  the  wall.  The  second  derivative  term 
is  important  only  in  the  region  near  the  wall.  Its  calculation 
is  quite  uncertain  and  is  given  mainly  to  show  its  order  of 
magnitude.  The  finite  intercept  at  the  wall  is  a  consequence 
of  the  arbitrary  linear  extrapolation  of  the  turbulence  inten¬ 
sity  to  zero.  'Pile  dissipation  term  is  the  most  diflieult  to 
determine  because  of  tin*  importance  of  the  region  near  the 
wall.  The  application  of  the  concept  of  local  isotropy  is 
strongly  conditioned  by  the  local  similarity  of  the  eddies 
responsible  for  the  dissipation,  the  verification  of  which  is 
dependent  on  the  dissipation  derivatives  obeying  the  iso¬ 
tropic  relation: 


y_ 

l 


Fioirk  2.T  Production  ami  dissipation  in  region  near  wall. 


Kioirk  2-1.  Production  and  dissipation  in  region  away  from  wall. 
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/&mV  .1 /jgV_  1  /d»V_ 

W/  2  2  vW _ 


From  figures  13  and  14  it  is  seen  that  only  in  the  outer 
region  of  the  boundary  layer,  f//4>0.7,  is  this  condition  satis¬ 
fied  and  that  the  derivatives  become  increasingly  divergent 

as  the  wall  is  approached.  The  terms  and  are 


equal  across  the  layer,  and  at  y/6~ 0.005  their  ratio  to  (£)' 

is  10  times  that  given  by  isotropy.  This  illustrates  the  very 
small  scale  nature  of  the  turbulent  motions  in  the  transverse 
directions  as  compared  with  the  longitudinal  in  the  region 
close  to  the  wall.  The  inadequacy  of  local  isotropy  is  strik¬ 
ingly  seen  from  the  overall  energy  balance  obtained  by  inte¬ 
grating  the  energy  equation  across  the  boundary  layer.  The 
integrals  of  the  diffusion  terms  and  the  second -derivative 
term  are  each  zero,  and  the  total  production  is  balanced  by 
the  total  dissipation  plus  the  total  convection.  Tsing  the 


isotropic  relation 


'•hit 


for  the  rate  of  dissipation  of 


turbulent  energy,  the  result  is  out  of  balance,  with  the  dissi¬ 
pation  being  loo  low  by  a  factor  of  2.1.  Hence,  in  order  to 
obtain  an  accurate  measure  of  t lit*  dissipation,  it  is  necessary 
to  measure  all  nine  derivatives.  Since  it  was  possible  to 
measure  only  the  first  five,  the  remaining  four  were  obtained 
by  the  arbitrary  assumption  that  derivatives  with  respect  to 
y  and  -  are  given  by  the  isotropic  relations 


/  dit  y  /  d r  y  /  dir  \a 
UyJ  -\6y)  \dy) 


/d«V  /dr  \‘  .j  dir  y 
(dr)  (dr)  \  dr  ) 


i  dr 


die 


'Phis  considerably  improved  the  total  energy  balance.  As 
seen  from  figure  23,  the  dissipation  calculated  in  this  manner 
was  obtained  as  close  to  the  wall  as  y; '5  =-0.005,  and  if  it  is  j 


Florae  26.  -Diffusion  and  convection. 


assumed  in  order  to  close  in  the  dissipation  curve  that  the 
dissipation  is  equal  to  the  production  in  the  region  (1<  y  4-. 
0.01,  the  disparity  is  15  percent.  This  ran  represent  a  sig¬ 
nificant  difference  in  the  )>oint-to-|M)int  balance  as  manifested 
by  the  diffusion  term.  When  the  diffusion  term  is  compared 
with  that  obtained  by  Townsend  (ref.  ti),  the  same  general 
trend  is  observed  although  there  is  a  considerable  difference 
in  magnitude.  The  necessary  balance  cannot  be  achieved 
by  the  gain  in  energy  due  to  diffusion  in  the  region  y  4  (Hi 
unless  it  is  cxtra|>nlatcd  to  an  unreasonably  long  distance. 
Consequently  the  conclusion  is  that  the  turbulent-energy 
dissipation  is  greater  than  the  production  in  a  very  thin 
region  next  to  the  wall  in  order  to  provide  the  dissipation 
needed.  This  means  a  diffusion  of  energy  toward  the  wall 
Townsend’s  values  for  the  kinetic-energy -diffusion  term 
require  a  movement  of  kenetic  energy  away  from  the  wall. 
Therefore,  if  energy  goes  toward  the  wall  it  must  be  due  to 
the  action  of  pressure  forces  (pressure  diffusion'.  This 
intense  dissipative  region  for  the  turbulence  apparently 
coincides  with  the  viscous  region  for  the  mean  flow.  The 
magnitude  of  the  diffusion  term  illustrates  the  weakness  of 
the  mixing-length  theory  which  has  the  implicit  assumption 
that  what  is  locally  produced  is  locally  dissipated. 

'Phis  flow  of  energy  against  the  energy  gradient  i-  rather 
surprising  and  contrary  to  what  is  intuitively  expected. 
It  does  not  seem  possible  that  such  a  result  could  be  brought 
about  incorrectly  by  the  questionable  assumption  made 
regarding  those  dissipation  derivatives  that  could  not  lie 
measured.  Neither  does  it  seem  reasonable  that  these 
derivatives  can  be  distributed  in  such  a  manner  as  to  make 
the  production  equal  to  the  dissipation  at  all  (Miints  across 
the  layer.  However,  it  is  felt  that  the  situation  is  compli¬ 
cated  by  the  phenomenon  of  intcrmitteney  and  it  would  lie 
helpful  to  have  a  fuller  understanding  of  the  role  that  it 
plays. 

The  energy  balance  for  the  mean  motion  is  conveniently 
expressed  in  terms  of  tin-  energy  thickness  4,  such  that 


where 


i v l' ‘X +" ('l s’  "i 


III 


Similarity  of  the  mean  velocity  distribution  permits  treating 
the  integral  in  the  expression  for  4,  as  a  constant.  If  4  is 
given  by 

4— ().37»l,4f 


then 


£-«wiv f.'r,  [■-('•-  )’]'%<» 


Equation  (1)  states  that  part  of  the  loss  of  kinetic  energy 
of  the  mean  motion  goes  directly  into  heat  through  the 
action  of  viscosity  and  the  remainder  goes  into  the  production 
of  turbulent  energy.  The  energy  balance  was  found  to  lx* 
satisfied  to  within  a  few  percent  and  served  as  a  welcome 
cheek  on  the  accuracy  of  the  measurements.  The  production 
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of  turbulent  energy  and  tin*  viscous  dissipation  arc  compared 
for  the  region  near  the  wall  in  figure  22.  This  emphasizes 
the  importance  of  the  region  near  the  wall.  In  fuel,  almost 
40  percent  of  the  loss  in  kinetic  energy  of  the  motion  is* 
directly  dissipated  by  viscosity  in  the  region  0<i/*<H0,  and 
of  the  remaining  tiO  percent  which  is  converted  into  turbulent 
energy  HO  percent  is  produced  in  the  same  region.  Thus,  if 
the  conclusion  drawn  from  the  turbulent-energy  balance  is 
correct.  So  percent  of  the  total  dissipation  (viscous  and 
turbulent)  takes  place  in  y*<H(). 

ENERGY  AND  SHEAR-STRESS  SPECTRA 

It  is  possible  to  gain  some  further  insight  into  shear  How 
by  examining  the  spectral  distributions  of  turbulent  energy 
and  shear  stress.  A  significant  advance  in  dealing  with  the 
energy  spectrum  has  been  made  in  the  domain  of  homo¬ 
geneous  and  isotropic  turbulence  (refs.  9  and  10).  The 
basic  concept  underlying  this  advance  is  that  energy  enters 
the  spectrum  through  the  large  eddies  and  is  then  trans¬ 
ferred  through  the  spectrum  to  the  smaller  eddies  where  it 
is  finally  dissipated.  If  the  lower  wave  numbers  are  excluded 
there  exists  a  range  in  which  the  eddies  are  in  a  state  of 
equilibrium,  governed  by  the  rate  at  which  they  transfer 
and  dissipate  energy.  When  the  Reynolds  number  is  high 
enough,  inertial  forces  will  predominate  in  the  lower  wave 
numbers  of  this  equilibrium  range,  and  a  relatively  pure 
transfer  region  will  exist.  By  dimensional  reasoning  it  can 
lie  shown  that  the  spectrum  will  vary  as  k~sn  in  this  range. 
Heisenberg  (ref.  11)  extended  this  concept  by  assuming 
that  the  transfer  of  energy  at  wave  number  k  was  caused 
by  a  turbulence  friction  produced  by  eddies  with  wave 
number  greater  than  /■.  He  represented  the  energy  balance 
in  the  equilibrium  range  for  homogeneous  and  isotropic 
turbulence  as 

iik"\  X*  2(m’<*')  '/*■'  (2) 

The  second  term  within  the  brackets  represents  a  turbulent 
viscosity  and  X  is  an  absolute  constant.  The  solution  of 
this  equation  gives  the  spectrum  in  the  following  form 

**  (:>  '•’  D<)7' 

where  k,  is  a  wave  number  in  the  intermediate  range  given  by 


For  low  wave  numbers  E(k)  varies  as  k~!,r\  and  for  the  high- 
wave-nttmber  end  where  viscous  forces  predominate,  as  k~‘. 

Although  equation  (2)  involves  an  assumed  mechanism 
that  may  not  entirely  represent  the  facts,  it  seems  to  be  a 
reasonably  good  approximation.  Where  the  Reynolds 
number  is  sufficiently  high,  there  is  evidence  of  a  J range 
and  a  transition  to  higher  negative  pow-ers  approaching  —7 
as  k  increases.  There  is  some  doubt  as  to  whether  —7  is 
the  correct  value  in  the  limit.  The  concept  embodied  in 
equation  (2)  provides  a  rational  basis  for  an  approach  to 


the  problem  of  shear  turbulence.  However,  the  extension 
to  shear  flow  is  complicated  by  the  presence  of  such  factors 
as  production,  diffusion,  convection,  and  the  absence  of 
isotropy  and  homogeneity.  Any  conclusions  as  to  the 
effects  of  the  diffusion  and  convection  are  difficult  to  draw. 
It  may  lie  assumed  that  such  effects  are  confined  to  the 
very  low  wave  numbers  which  lie  outside  the  equilibrium 
range.  In  addition,  the  convection  may  be  considered 
negligible  across  most  of  the  layer.  An  attempt  to  assess 
thi‘  influence  of  the  production  term  in  the  equilibrium 
range  of  the  spectrum  was  made  by  Tchcn  (ref.  2H).  By 
considering  the  influence  of  the  mean  velocity  gradient  the 
conclusion  is  reached  that  a  range  of  k~'  will  exist  in  the 
wave-number  region  where  k~:''  normally  exists  when  there 
is  no  gradient. 

An  experimental  tost  of  theoretical  predictions  is  rendered 
difficult  by  the  fact  that  only  the  one-diimuisiona)  spectrum 
can  be  measured.  In  isotropic  turbulence  this  is  not  too 
restrictive  a  factor  because  the  transformation  from  the 
three-dimensional  spectral  function  such  as  appears  in 
equation  (H)  to  the  one-dimensional  spectral  function  is 
known  (ref.  11).  In  shear  flow  this  relation  is  not  known, 
and  one  has  to  lx*  content  with  the  qualitative  inference  that 
in  some  unknown  manner  the  one-dimensional  spectrum  is 
still  an  integral  effect  of  the  three-dimensional.  Despite  the 
aforementioned  complications  tin1  measured  spectra  arc  of 
interest,  and  several  interesting  features  can  be  noted. 

The  spectra  of  u!  at  various  cross-sectional  positions  are 
given  in  figures  <i  and  7.  The  trend,  in  going  toward  the 
surface,  is  for  the  higher  wave  numbers  to  have  a  greater 
percentage  of  the  turbulent  energy.  This  is  in  accord  with 
the  trend  of  the  shear-stress  spectra  (fig.  X)  also  to  extend  to 
higher  and  higher  wave  numbers  as  the  wall  is  approached. 
It  is  noted  that  nearly  all  of  the  turbulent  energy  lies  within 
the  stress-producing  range.  The  spectrum  of  shear  stress  at 
yl&—  0.2  shows  an  increase  over  that  at  y  4=0.00  in  the  lower 
wave  numbers.  This  may  be  indicative  of  the  influence  of 
the  wall  becoming  negligible  at  t/6 -0.2.  For  y  6)().J  the 
decrease  in  shear  stress  across  the  boundary  layer  takes  place 
for  the  entire  spectral  range.  All  of  tin1  energy  spectra 
indicate  the  existence  of  a  region  varying  as  k\~~  at  the  high- 
number  end.  but  because  of  the  lesser  accuracy  of  the  meas¬ 
urements  in  this  range  no  direct  comparison  has  been  made. 
In  the  outer  region  of  the  layer  at  y  6-  0.58  and  0.8  there  is 
an  extensive  region  where  the  spectrum  of  ir  varies  as  ;l 
corresponding  to  the  inertial  subrange.  Tin*  effect  of  inter- 
rnitteney  in  this  region  is  difficult  to  ascertain  except  to  say 
that  it  may  be  confined  to  the  low-wave-number  end  of  the 
spectrum.  Also  the  maximum  in  the  spectrum  at  »/  4-  1.0  is 
apparently  a  consequence  of  intermit tency.  There  is  a 
gradual  transition  in  the  shape  of  the  spectrum  front ;/  6-0.58 
to  t//6=0.05  where  there  is  a  wave-number  range  with  the 
slope  kr‘  as  predicted  by  Tclien.  The  small-scale  nature'  of 
the  turbulence'  near  the  wall  is  shown  by  the  spectrum  of  uJ 
at  y/6=  0.00 1 1 .  The1  elip  in  the  spectrum  at  the  low-wave- 
numltcr  erne!  as  inelie-ntcel  by  the  elasheel  curve  may  l>c  elite  to 
e'.xpe'ri mental  rrrttr.  Howe'ver.  it  was  repeatable  anel  at  ten- 
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tiim  is  drawn  to  it  because  it  may  la*  a  result  of  some  charae- 
teristie  phenomenon  associated  with  the  laminar  sublayer. 

Empirically,  the  isotropic  relation  for  the  dissipation  is  in 
fair  agreement  with  that  given  in  figure  24  from  >/  &  0.05  to 

1.0.  fly  assuming  that  the  second  moments  of  the  spectra 
(fig.  20)  are  a  fair  representation  of  the  total  dissipation,  it 
is  seen  that  in  going  toward  the  wall  the  shear  penetrates 
deeper  and  deeper  into  the  dissipation  spectrum.  For 
example,  at  ;/ '{  —  0.58  the  wave-number  range  up  to  where 


L  has  fallen  to  10  *  contains  approximately  It)  per- 

■  i 

cent  of  the  dissipation.  At  y,  5 -0.05  this  figure  is  TO  percent. 
It  is  reasonable  to  export  this  trend  to  continue,  with  the 
r<‘sult  that  very  close  to  the  wall  the  turbulence  produced  is 
directly  dissipated  and  there  is  no  significant  transfer  of 
energy  along  the  spectrum.  This  conclusion  is  consistent 
with  the  small-scale  turbulence  near  the  wall  and  the  rapidly 
decreasing  spectrum  at  ij  5—0.0011. 

The  shear  spectrum  is  a  direct  test  of  local  isotropy,  and 
it  is  evident  that  the  range  of  wave  numbers  for  which  local 
isotropy  exists  becomes  progressively  smaller  as  the  wall  is 
approached.  The  transfer  range,  as  evidenced  by  t',  ,v\  bus 
not  yet  become  locally  isotropic  since  the  shear  spectra 
approach  zero  at  some  higher  wave  number.  However,  the 
energy  spectrum  is  apparently  insensitive  to  the  small  amount 
of  shear  stress  which  does  exist.  In  figures  !)  and  10  the 
measured  spectra  of  iJ  at  y/6  0.05  and  0.58  are  compared 

with  those  calculated  from  the  measured  it-  spectra  using  the 
isotropic  relation 


,1 

2  r Ik, 


/»<*,) 


which,  strictly  speaking,  is  valid  only  for  the  range  of  local 
isotropy.  It  is  seen  that  the  measured  r-  spectra  agree  with 
the  calculated  at  the  higher  wave  numbers.  This  is  con¬ 
sistent  with  the  range  of  local  isotropy  indicated  b\  the 
spectra  of  shear,  and  the  difference  between  the  overall  <i’ 
and  r  is  confined  to  the  lower  wave  numbers. 

From  the  turbulent-energy  equations  for  the  individual 
velocity  components  it  is  seen  that  the  terms  responsible  for 
the  exchange  of  energy  among  the  various  components  are 


.  du  dr 

the  pressure  terms  p,  p ,  »  and  />, 

of  such  transfer  is  unknown  but  a 


■  The  speet  ral  range 
reasonable  speculation 


would  be  that  the  pressure  term  p,  ^  responsible  for  the 

transfer  of  energy  to  r-  exists  in  the  same  wave-number 
range  as  the  shear-stress  spectra.  There  is  little  evidence 
in  the.  measured  r-  spectrum  of  a  range  of  k, " 1  as  observed 
for  the  spectrum  of  This  may  be  a  consequence  of  the 
spectral  distribution  of  the  pressure  term  being  different 
from  that  of  the  shear  stress  with  a  weighting  more  to  the 
higher  wave  number  end  of  ibis  rang*'. 


CONCLUSIONS 

Tin-  following  conclusions  were  obtained  from  an  investi¬ 
gation  of  charnel  eristics  of  turbulence  in  a  boundary  layer 
with  zero  pressure  gradient.  Here  y  is  the  distance  normal 
to  the  surface  measured  from  the  surface.  »/*  is  y  times  the 
square  root  of  the  shearing  stress  at  the  wall  divided  by  tin1 
density  of  air  over  the  kinematic  viscosity  i/f\V  (where 
I',  \  ro/p),  and  6  is  the  boundary-layer  thickness. 
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1.  The  turbulent  boundary  layer  ran  lie  considered  to  be 
divided  into  three  regions:  A  viscous  region  extending 
beyond  the  laminar  sublayer  to  a  value  of  //*  —  .’{(),  an  inter¬ 
mediate  region  from  ;/*  =  :{()  to  approximately  ;/  <5  =  0.2  where 
the  influence  of  the  wall  still  exists,  and  an  outer  region 
characterized  by  the  phenomenon  of  intermit tenev. 

2.  The  turbulent  boundary  layer  has  a  sharp  outline 
constantly  changing  in  an  irregular  manner.  The  turbulent  j 
regions  travel  with  a  velocity  lower  than  that  of  the  free  j 
stream.  A  rough  estimate  of  the  average  wave  length  of  this  i 
irregular  boundary  would  he  approximately  2 6. 

d.  The  position  of  the  edge  of  the  layer  with  respect  to  I 
the  surface  is  approximated  closely  by  a  Gaussian  distribution 
centered  at  0.7X6  with  a  standard  deviation  equal  to  0.146. 

4.  The  probability  density  of  the  a-lltictuation  is  very 
close  to  Gaussian  in  the  region  >/  6<4).4. 

">.  The  importance  of  the  region  near  the  wall  has  been 
demonstrated;  in  fact,  approximately  X*>  percent  of  the  total 
dissipation  (viscous and  turbulent )  takes  place  within 
The  production  of  turbulent  energy  and  the  turbulent -energy 
dissipation  reach  a  sharp  maximum  within  this  region,  and 
there  is  an  inward  How  of  energy  toward  the  wall  because  of 
the  action  of  pressure  forces. 

f>.  The  concept  of  local  isotropy  is  inadequate  for  obtaining  1 
the  turbulent -energy  dissipation,  cspeciallv  in  the  region  near  ! 
the  wall. 

7.  In  the  outer  region  of  the  boundary  layer  where  the 
mean  velocity  gradient  is  small  the  existing  theories  for  the 
energy  spectrum  in  isotropic  turbulence  can  be  applied. 
As  the  wall  is  approached  the  nonisotropy  becomes  significant 
with  a  resultant  change  in  the  spectrum.  Very  close  to  the 
wall,  the  turbulent  energy  produced  is  directly  dissipate!  with 
no  significant  transfer  of  energy  through  tin*  spectrum. 

X.  As  dose  to  the  wall  as  »/  6  0.0.')  it  has  been  shown  that 

the  space-time  transformation  is  valid  for  the  smaller  eddies 
responsible  for  the  turbulent  dissipation  but  becomes  invalid 
for  tin*  larger  eddies. 

National  Hi  heap  ok  Standards. 

Washington,  I).  (’.,  Min/ S,  I !>')■{. 

REFERENCES 

1.  Corrsin.  Stanley:  Invest  illation  of  How  in  an  Axially  Svitnet- 

rical  Heated  Jet  of  Air.  \ACA  W|{  W  !M.  1943.  (Formerly 

NACA  AC  It  31,2 3J 

2.  'fownseiid,  A.  A.:  The  Fullj  Developed  Turbulent  Wake  of  a 

Circular  Cylinder.  Australian  Jour.  Sci.  He?*..  ser.  A,  vol.  2.  ! 

no.  I.  Dee.  1949,  pp.  451  Ills. 

3.  Townsend.  A.  A.:  Local  Isotropy  in  the  Turbulent  Wake  of  a  ! 

Cylinder.  Australian  Jour.  Sri.  He*..  ser.  A.  vol.  I,  no.  2. 

June  1948,  pp.  Iff]  1 7 1. 

4.  Laufer.  John:  Investigation  of  Turbulent  Flow  in  a  Two-Dimen¬ 

sional  Channel.  NACA  Hep.  1 053.  1951 .  (Supersedes  XACA 

TN  2123.) 


a.  Laufer.  John:  The  Strueture  of  Turbulence  in  Fully  Developed 
Pipe  Flow.  NACA  Hep.  117-4.  1954.  (Supersede-*  NACA 
TX  21154. i 

0.  Townsend,  A.  A.:  The  Strueture  of  the  Turbuleni  lioundarx 
layer.  Proe.  Cambridge  Phil.  Sor.,  vol.  17.  pt.  2.  Apr.  1951. 
pp.  375  395. 

7.  Sehubauer.  (I.  H.,  and  Klebanoff.  P.  S.:  Investigation  « *f  Separa¬ 
tion  of  the  Turbulent  Houndary  Layer.  NACA  Hep.  1030, 
1951.  (SujMTsedes  XACA  TX  2133.) 

K.  KlebanolT.  P.  S..  and  Diehl.  Z.  W.:  Some  Features  of  Artificially 
Thickened  Fully  Developed  Turbulent  Boundary  Layers  With 
Zero  Pressure  (Iradienl..  XACA  Hep.  II  It).  1952.  iSuper- 
sedes  X  AC  A  TX  2  475.) 

9.  Kolmogoroff.  A.  X.:  The  Loeal  Strueture  of  Turbuienee  in  Incom¬ 
pressible  Fluid  for  \‘ery  Large  Reynolds'  Number-.  Comp 
rend..  Acad.  Sci.  CUSS,  vol.  30.  no.  4.  Feb.  10,  1941,  pp.  301  305. 

10.  Kolmogorotr,  A.  X.:  Dissipation  of  Knergy  in  the  Locally  Iso¬ 

tropic  'rurbulenee.  Comp,  rend.,  Aead.  Sci.  CHSS.  vol.  32. 
no.  1,  July  10,  1911,  pp.  Hi  IS. 

11.  Heisenberg.  W\:  Zur  st  at  ist  ischen  Theorie  der  Turbulenz.  Z-. 

Phvs..  Bd.  124.  Heft  7  12.  19  IS.  pp.  028  #>57. 

12.  Kovasznay.  Leslie  S.  <«.:  Development  of  Turbulence-Measuring 

equipment .  XACA  Hep.  1209.  1954.  ‘Supersede-*  XACA 

TX  2839. ) 

13.  Sehubauer.  (I.  B.,  and  KlebanolT.  P.  S.:  Theory  and  Application 

of  Hot-Wire  Instruments  in  the  Investigation  of  Turbulent 
Boundary  Layers.  XACA  WH  W  Sti.  1940.  i  Formerly 
XACA  AC  It  5K27.) 

II.  Xewman,  B.  (1..  and  Leary.  B.  O,:  The  Measurement  of  tin* 
Reynolds  Stresses  in  a  Circular  Pipe  as  a  Means  of  Testing  a 
Hot  Wire  Anemometer.  Hep.  A. 72,  Dept,  of  Supply.  Aero 
Res.  Lab..  Australia.  Xov.  1950. 

15.  Squire.  II.  B..  and  Young.  A.  I).:  The  Caleulat ion  of  the  Profile 
Drag  of  Aerofoils.  R.  A*  M.  Xo.  IS3S.  British  A.  R.  C.,  1938. 
Hi.  Liepmann.  II.  W  ..  and  Robinson.  M.  S.:  ('minting  Methods  and 
Fquipment  for  Mean-Value  Measurements  in  Turbulence 
Research.  XACA  TX  3037.  1953. 

17.  'Taylor.  (1.  L:  Statistical  Theory  of  Turbulence  III  Distribution 
of  Dissipation  of  Knergy  in  a  Pipe  Over  Its  Cross-Section. 
Proc.  Roy.  Sor.  (London),  ser.  A,  vol.  151,  no.  873.  Sept.  2. 
1935,  pp.  455  404. 

IS.  Dryden.  Hugh  I...  Sehubauer.  (i.  B..  Mock.  W.  C..  Jr.,  and 
Skrairstad.  H.  K.:  Measurements  of  Intensity  and  Scale  of 
W uid-Tuunrl  Turbulence  ami  Their  Relation  to  the  Critical 
Reynolds  Nun  ber  of  Spheres.  XACA  Rep.  5SI.  1937. 

19.  Cberoi.  Mahinder  S..  and  Kovasziiay.  Leslie  S.  (1.:  Influence  of 

Resolving  Power  on  Measurement  of  Correlations  ami  Spectra 
of  Random  Fields.  Tech.  Hep.  No.  30.  Project  Squid.  John- 
Hopkins  Cniv.,  1952. 

20.  Kovasznay,  L.  S.  (’».:  Quarterly  Progress  Report  of  Aeronautics 

Department.  Contract  Nord  S030  JHB  31).  Bur.  Ord.,  Dept 
Navy,  and  'Die  Johns  Hopkins  Cniv..  Jan.  I  Mar.  31.  1950. 

21.  Corrsin,  Stanley,  and  Kistler.  Alan  !..:  'The  Free-Slream  Bound¬ 

aries  of  Turbulent  Flows.  XACA  Hep.  1244.  1955.  (Super¬ 
sedes  XACA  TX  3133.) 

22.  (Virrsin,  S.:  Interpretation  of  Viscous  'Terns  in  the  Turbulent 

Knergy  Kquntinn.  Readers’  Forum.  Jour.  Aero.  Sci..  vol.  20, 
no.  12.  Dec.  1953.  pp.  S53  854. 

23.  Tchen.  C.  M.:  On  the  Spectrum  of  Knergy  in  'Turbulent  Shear 

Flow.  Res.  Pu|M»r  10*2388,  Jour.  Res.,  Xat.  Bur.  Standards, 
vol.  50.  no.  I.  Jan.  1953.  pp.  51  02. 


u.  <l  60¥I«R«fNT  FfttHTINE  OrriCT  ttlft 


